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I. INTRODUCTION 



The two-photon coupling of axions or axion-like particles allows for transitions between them and photons in the 
presence of external electric or magnetic fields as shown in Figure 1 0, 01 • This mechanism serves as the basis for the 
experimental searches for galactic dark matter axions 1.,..3] and solar axions EL 0, IE IS • The astrophysical implica- 
tions of this mechanism have also been widely investigated and reviewed 0, Q • The phenomenological consequences 
of an extremely light or massless axion would be particularly interesting in several astrophysical circumstances such 
as polarization of radio-galaxies |lf]j | and quasars the diffuse x-ray background |T^ . or ultra-high energy cosmic 
rays [3[l|. 

One intriguing cosmological consequence of this mechanism is photon-axion conversion caused by intergalactic 
magnetic fields, leading to the dimming of distant light sources, notably of supernovae (SNe) of type la that are used 
as cosmic standard candles Observationally, SNe la at redshifts 0.3 < z < 1.7 appear fainter than expected from 
the luminosity-redshift relation in a decelerating Universe [Hi Il7l Il8| , a finding usually interpreted as evidence for 
acceleration of the cosmic expansion rate and thus for a cosmic equation of state (EoS) that today is dominated by 
a cosmological constant, a slowly evolving scalar field, or something yet more exotic [l9j . The dimming caused by 
photon-axion conversion could mimic this behavior and thus provide an alternative to the interpretation as cosmic 
acceleration. Although still requiring some non-standard fluid to fit the flatness of the universe, this model seemed 
capable to explain the SN dimming through a completely different mechanism. 

However, if the light from distant SNe la reaches us partially converted to axion-like particles, the same mechanism 
would affect all distant sources of electromagnetic radiation. Therefore, it appears useful to update the different 
arguments constraining photon-axion conversion in intergalactic magnetic fields, in particular the constraints arising 
from spectral distortions of the cosmic microwave background (CMB) and dispersion of quasar (QSO) spectra. 

To this end we begin in Section 2 with a review of the formalism of photon-axion conversion in magnetic fields. 
Some technical details are deferred to Appendix A. In Section 3 we describe how this mechanism affects the SN la 
luminosity-redshift relation and accounts for the observed dimming. In Section 4 we turn to spectral CMB distortions 
and in Section IV! combine these limits with those from the dispersion of QSO spectra. In Section 6 we describe some 
additional limits from a violation of the reciprocity relation between the luminosity and angular diameter distances. 
We conclude in Section IVIII and comment on the viability of the photon-axion conversion mechanism. 



II. PHOTON-AXION CONVERSION 



To understand how photon-axion conversion could affect distant sources we take a closer look at the phenomenon 
of photon-axion mixing. The Lagrangian describing the photon-axion system is |8j 



~\~ C*a "F 



(1) 



The QED Lagrangian for photons is 



r — F F^ v 



90mi L 



(2) 



F pa its dual, a the fine-structure constant, and 



where F^v is the electromagnetic field-strength tensor, F^ v — ^e^ vpa 
m e the electron mass. We always use natural units with h = c = kn = 1. The second term on the r.h.s. is the 
Euler-Heisenberg effective Lagrangian, describing the one-loop corrections to classical electrodynamics for photon 
frequencies u> <C m e . 




FIG. 1: Axion-photon transition in an external electric or magnetic field. 
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The Lagrangian for the non-interacting axion field a is 

£a = -d^ad^a - ^m 2 a 2 . (3) 

A generic feature of axion models is the CP-conserving two-photon coupling, so that the axion-photon interaction is 

£ ai = -\ga-,F^F^a = g ai E B a, (4) 

where g ai is the axion-photon coupling with dimension of inverse energy. A crucial consequence of C is that the 
propagation eigenstates of the photon-axion system differ from the corresponding interaction eigenstates. Hence 
interconversion takes place, much in the same way as it happens for massive neutrinos of different flavors. However, 
since the mixing term F^F^a involves two photons, one of them must correspond to an external field pi 13. I§1 120| . 

Axion-photon oscillations are described by the coupled Klein-Gordon and Maxwell equations implied by these 
Lagrangians. For very relativistic axions [m a <C to), the short- wavelength approximation can be applied and the 
equations of motion reduce to a first order propagation equation. More specifically, we consider a monochromatic light 
beam travelling along the z-direction in the presence of an arbitrary magnetic field B. Accordingly, the propagation 
equation takes the form 



{u-id z +M)\ A y =0, (5) 




where A x and A y correspond to the two linear photon polarization states, and u) is the photon or axion energy. The 
mixing matrix is 

(^XX A x y "n,QayB x \ 

^y* ^yy b«~< B y . ( 6 ) 

where A a = —m 2 1 /2uj. The component of B parallel to the direction of motion does not induce photon-axion mixing. 
While the terms appearing in the third row and column of A4 have an evident physical meaning, the Ay -terms 
(i,j = x,y) require some explanations. Generally speaking, they are determined both by the properties of the 
medium and the QED vacuum polarization effect. We ignore the latter, being sub-dominant for the problem at 
hand |21|. 

For a homogeneous magnetic field we may choose the y-axis along the projection of B perpendicular to the z-axis. 
Correspondingly we have B x = 0, B y — |Bt| = B sm8, A x — A±, A y — A\\. Equation © then becomes 



(uj-id z +M) ( A j ) =o. (7) 

with the mixing matrix 






A R 









A a7 





A a7 


A a 



M = | \,< A A„- 1 . (M 

where 



A Q7 = g a7 |B T |/2 , (9) 
V 



A, u = A pl + A™ . (10) 



In a plasma, the photons acquire an effective mass given by the plasma frequency uj 2 ^ = Aitan e /m e with n e the 
electron density, leading to 

Furthermore, the A™ terms describe the Cotton-Mouton effect, i.e. the birefringence of fluids in the presence of a 
transverse magnetic field where |A|? M — A^ M | oc B\. These terms are of little importance for the following arguments 
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and will thus be neglected. Finally, the Faraday rotation term Ar, which depends on the energy and the longitudinal 
component B z , couples the modes and Aj_. While Faraday rotation is important when analyzing polarized sources 
of photons, it plays no role for the problem at hand. 

With this simplification the Aj_ component decouples, and the propagation equations reduce to a 2-dimensional 
mixing problem with a purely transverse field B = Bt 



(u J -id z +M 2 ) ( A|1 ) =0, (12) 



with a 2-dimcnsional mixing matrix 



The solution follows from diagonalization by the rotation angle 

In analogy to the neutrino case |22j . the probability for a photon emitted in the state A\\ to convert into an axion 
after travelling a distance s is 

i%( 7 -o) = |(A||(0)|a( S )}| 2 

= sin 2 (2?9)sin 2 (A osc s/2) 
2 sin 2 (A osc s/2) 

= (A - S) (A oscS /2) 2 < (15) 



where the oscillation wavenumber is given by 

A 2 sc = (A pl -A a ) 2 +4A 2 7 . (16) 

The conversion probability is energy-independent when 2|A a7 | ^ |A p i — A a | or whenever the oscillatory term in 
Eq. Ijl5|l is small, i.e. A OS cs/2 <C 1, implying the limiting behavior Pq = (A Q7 s) 2 . 

The propagation over many i?-field domains is a truly 3-dimensional problem, because different photon polarization 
states play the role of An and A± in different domains. This average is enough to guarantee that the conversion 
probability over many domains is an incoherent average over magnetic field configurations and photon polarization 
states. The probability after travelling over a distance r 3> s, where s is the domain size, is derived in Appendix A 
along the lines of Ref. |2jj and is found to be 



1 — exp 



3P r 
2s 



(17) 



with Pq given by Eq. (|15|) . As expected one finds that for r/s — > oo the conversion probability saturates, so that on 
average one third of all photons converts to axions. 

III. PHOTON-AXION CONVERSION AND SUPERNOVA DIMMING 

A. Observations 

In 1998, two groups using SNe la as cosmic standard candles reported first evidence for a luminosity-redshift relation 
that indicated that the expansion of the universe was accelerating today [161 Il7| . The quantity relevant for SN la 
observations is the luminosity distance 4 at redshift z, defined by 

dj (z) = — , (18) 

where C is the absolute luminosity of the source and T is the energy flux arriving at Earth |l6l Il7| . In Friedmann- 
Robertson- Walker cosmologies, the luminosity distance at a given redshift z is a function of the Hubble parameter 
.Ho, the matter density Om, and the dark energy density Oa- Usually the data are expressed in terms of magnitudes 

m = M + 51og 10 (J±^j +25, (19) 
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FIG. 2: SN la Hubble diagram. Upper panel: Hubble diagram for low- and high-redshift SN la samples. Overplotted are three 
cosmologies: "low" and "high" S!m with = and the best fit for a flat cosmology Cm = 0.24 and Q\ = 0.76. Lower panel: Difference 
between data and models with Qm = 0.20 and Qa = 0. (Figure from Ref. |lrJl with permission.) 



where M is the absolute magnitude, equal to the value that m would have at = 10 pc. 

Figure shows the Hubble diagram for SN la samples at low and high z. The distances of high-redshift SNe are, 
on average, 10% to 15% larger than in a low matter density (Om = 0.2) Universe without dark energy (Oa = 0). 
Therefore, objects of a fixed intrinsic brightness appear fainter if the cosmic energy density budget is dominated by 
dark energy. The best fit of these data supports a Universe composed of a fraction of dark matter Qm — 0.3 and a 
fraction of dark energy fl\ ~ 0.7. 

Dark energy has been associated with vacuum energy or an Einstein cosmological constant that produces a constant 
energy density at all times. Defining the equation of state 

w = V - , (20) 
P 

the cosmological constant is characterized by p = —p, i.e. w = — 1. From the Friedmann equations any component of 
the density budget with equation of state w < —1/3 causes cosmic acceleration. SN la data imply that w > —0.5 are 
disfavoured, supporting the cosmic acceleration of the Universe |Tt| . 



B. Interpretation in terms of photon-axion conversion 



To explore the effect of photon-axion conversion on SN dimming we recast the relevant physical quantities in terms 
of natural parameters. The energy of optical photons is a few eV. The strength of widespread, all-pervading S-fields 
in the intergalactic medium must be less than a few 10~ 9 G over coherence lengths s crudely at the Mpc scale, 
according to the constraint from the Faraday effect of distant radio sources [24|. Along a given line of sight, the 
number of such domains in our Hubble radius is about N » H^ 1 / s w 4 x 10 3 for s ~ 1 Mpc. The mean diffuse 
intergalactic plasma density is bounded by n e < 2.7 x 10 -7 cm -3 , corresponding to the recent WMAP measurement of 
the baryon density pE| . Recent results from the CAST experiment give a direct experimental bound on the axion- 
photon coupling of g ai ^ 1.16 x 10~ 10 GeV -1 , comparable to the long-standing globular-cluster limit 8J. Suitable 
representations of the mixing parameters are 

0.15 5io B nG , 

-""""(^"(tS?)" 1 . 
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FIG. 3: Hubble diagram for SNe la for different cosmological models, relative to the curve with Otot = (dotted horizontal line). The 
dashed curve is a best fit to the SN data assuming that the Universe is accelerating (S!m = 0.3, Qjy = 0.7); the solid line is the photon-axion 
oscillation model with fl^ = 0.3 and Q$ = 0.7, the dot-dashed line is SIm = 0.3, Qg = 0.7 with no oscillation, the dot-dot-dashed line is 
for Qj£ = 1 and again no oscillation. (Figure from Ref. Il r l with permission.) 



Mpc" 



= -11.1 



leV 



10 7 cm -3 



(21) 



where we have introduced gio — <? O7 /10 10 GeV 1 and B n Q is the magnetic field strength in nano-Gauss. 

The mixing angle defined in Eq. Ijl4(l is too small to yield a significant conversion effect for the allowed range of 



axion masses because IAJ » |A„ 



VI 



Therefore, to ensure a sufficiently large mixing angle one has to require 



nearly massless pseudoscalars, sometimes referred to as "arions" |26|,|27]]. For such ultra-light axions a stringent limit 
from the absence of 7 -rays from SN 1987A gives g ai < 1 x 1CT 11 GeV" 1 [H or even g aj < 3 x KT 12 GeV" 1 HU. 
Henceforth we will consider the pseudoscalars to be effectively massless, so that our remaining independent parameters 
are <7io-B n G and n e . Note that m a only enters the equations via the term m 2 — uj^, so that for tiny but non-vanishing 

values of m a , the electron density should be interpreted as rt e , c fr = \n e — rn^m e / '(47ra)|. 

Allowing for the possibility of photon-axion oscillations in intergalactic magnetic fields, the number of photons 



emitted by the source and thus the flux T is reduced to the fraction P 7 ^ 7 
distance [Eq. l(T%|l ] becomes 



d L 



d L /P^ 



1 — P 7 ^ a - Therefore, the luminosity 



(22) 



and the brightness [Eq. (|19fl ] 



™ ~ 2 lo glo( P 7^7)- 



(23) 



Distant SNe la would eventually saturate (P 7 ^ 7 = 2/3), and hence they would appear (3/2) x / 2 times farther away 
than they really are. This corresponds to a maximum dimming of approximately 0.4 mag. Csaki, Kaloper and 
Terning (CKT I) showed that if photon-axion conversion takes place, this mechanism can reproduce the SN Hubble 
diagram [lf| , assuming, for example, a nonstandard dark energy component Q$ — 0-7 with equation of state w = —1/3, 
which does not produce cosmic acceleration (Fig.|3J|. 

However, in the model of CKT I, plasma density effects were neglected (n e = 0). Later it was recognized that the 
conclusions of CKT I can be significantly modified when the effects of the intergalactic plasma on the photon-axion 
oscillations are taken into account [2l|. In the presence of plasma effects, the probability of oscillation is lower than 
before and it is no longer achromatic (Fig. SN observations not only require dimming, but also that the dimming 
is achromatic. In fact, SN observations put a constraint on the color excess between the B and V bands, 



E[B -V] = -2.5 log 



10 



F°{B) F e {V) 
F C {B) F°(V) 



(24) 



where F° or F e is the observed or emitted flux, respectively. The B and V band correspond to 0.44 /jm and 0.55 /im, 
respectively. Observations constrain E[B — V] to be lower than 0.03 01- This can be translated to 



P( 7 -> a) v 



Pin a)v 



- 1 



< 0.03 



(25) 
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FIG. 4: Ratio of the probability of conversion of photons to axions including the effects of the intergalactic plasma (n e ss 10 _7 cm — 3 ) 
and the probability of oscillations when this effect is not considered, as a function of the photon energy u). The curves are d raw n for 
different size s of the magnetic domains: 0.5 Mpc (dashed line), 1 Mpc (solid line) and 2 Mpc (dotted line). (Figure from Ref. L'll with 
permission.) 

Therefore, assuming an electron density n e ribaryons = n^rj ~ 10~ 7 cm~ 3 , the model is ruled out in most of the 
parameter space because of either an excessive photon conversion or a chromaticity of the dimming plj . Only fine- 
tuned parameters for the statistical properties of the extragalactic magnetic fields would still allow this explanation. 

On the other hand, Csaki, Kaloper and Terning [3(j (CKT II) criticized the assumed value of n e as being far too 
large for most of the intergalactic space, invoking observational hints for a value at least one order of magnitude 
smaller. As a consequence, for values of w p i < 6 x 10~ 15 eV, corresponding to n e < 2.5 x 10~ 8 cm~ 3 , one finds 
\Pv — Pb\ < 0.03 so that the chromaticity effect disappears very rapidly, and becomes undetectable by present 
observations. 

Figure shows qualitatively the regions of n e and gioB n Q relevant for SN dimming at cosmological distances. To 
this end we show iso-dimming contours obtained from Eq. (|23[1 for a photon energy 4.0 eV and a magnetic domain 
size s = 1 Mpc. For simplicity we neglect the redshift evolution of the intergalactic magnetic field B, domain size 
s, plasma density n e , and photon frequency to. Our iso-dimming curves are intended to illustrate the regions where 
the photon-axion conversion could be relevant. In reality, the dimming should be a more complicated function since 
the intergalactic medium is expected to be very irregular: there could be voids of low n e density, but there will also 
be high density clumps, sheets and filaments and these will typically have higher B fields as well. However, the 
simplifications used here are consistent with the ones adopted in CKT II and do not alter our main results. 

The iso-dimming contours are horizontal in the low-n e and low-f/io-BnG region. They are horizontal for any gioB n Q 
when n e is sufficiently low. From the discussion in Sec. [H] we know that the single-domain probability Pq of Eq. i|15|) 
is indeed energy independent when |A osc s| <§; 1, i.e. for |A p i|s/2 <C 1 and |A a7 |s <C 1. When n e < few 10 -8 cm -3 
and <?io-BnG ^5 4, we do not expect an oscillatory behavior of the probability. This feature is nicely reproduced by 
our iso-dimming contours. From Fig. [S] we also deduce that a significant amount of dimming is possible only for 
5io-BnG > 4 x 10" 2 . 

In CKT I, where the effect of n e was neglected, a value m a ~ 10~ 16 eV was used. In terms of our variables, 
this corresponds to n e , e ff « 6 x 10 -12 cm -3 . As noted in CKT II, when plasma effects are taken into account, any 
value n e < 2.5 x 10~ 8 cm" 3 guarantees the required achromaticity of the dimming below the 3% level between the 
B and V bands. The choice B n Q of a few and 1710 ~ 0.1 in CKT I and II falls in the region where the observed SN 
dimming could be explained while being marginally compatible with the bounds on the intergalactic B field and on 
the axion- photon coupling gig. 

IV. CMB CONSTRAINTS 

If photon-axion conversion over cosmological distances is responsible for the SN la dimming, the same phenomenon 
should also leave an imprint in the CMB. A similar argument was previously considered for photon- graviton con- 
version |3lj |. Qualitatively, in the energy-dependent region of P 7 ^ Q one expects a rather small effect due to the low 
energy of CMB photons (u> ~ 10~ 4 eV). However, when accounting for the incoherent integration over many domains 
crossed by the photon, appreciable spectral distortions may arise in view of the accuracy of the CMB data at the 
level of one part in 10 4 -10 5 . For the same reason, in the energy- independent region, at much lower values of n e than 
for the SNe la, the constraints on gioB n Q are expected to be quite severe. The depletion of CMB photons in the 
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FIG. 5: Iso-dimming curves for an attenuation of 0.01, 0.1, and 0.4 magnitudes. The photon energy of 4.0 eV is representative of the 
B-band. The size of a magnetic domain is s = 1 Mpc. (Figure from Ref. 1313 . ) 



patchy magnetic sky and its effect on the CMB anisotropy pattern have been previously considered [15j. However, 
more stringent limits come from the distortion of the overall blackbody spectrum [3^ | . 

To this end the COBE/FIRAS data for the experimentally measured spectrum were used, corrected for fore- 
grounds ■ Note that the new calibration of FIRAS [2] is within the old errors and would not change any of our 
conclusions. The N — 43 data points $° xp at different energies u>i are obtained by summing the best-fit blackbody 
spectrum (Fig. [fJJ) to the residuals reported in Ref. (3j| • The experimental errors Oi and the correlation indices 
between different energies are also available. In the presence of photon-axion conversion, the original intensity of the 
"theoretical blackbody" at temperature T 



^VTHf^exp^/T)-!]- 1 
would convert to a deformed spectrum that is given by 



In Ref. [32| , we build the reduced chi-squared function 



xl(T,X) 



N 



1 N 



(26) 



(27) 



(28) 
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FIG. 6: Uniform CMB spectrum and fit to the blackbody spectrum. Uncertainties are a small fraction of the line thickness. (Figure 
from Ref. \ '> j with permission.) 
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FIG. 7: Exclusion plot for axion-photon conversion based on the COBE/FIRAS CMB spectral data. The region above the solid curve 
is excluded at 95% CL whereas the one above the dotted curve is excluded at 99% C L. T he size of each magnetic domain is fixed at 
s = 1 Mpc. We also reproduce the iso-dimming contours from Fig.|3 (Figure from Ref. lolj .i 

where 

A$, = - *V,, T)P^{u u A) (29) 

is the i-th residual, and 

a % = PijViVj (30) 

is the covariance matrix. We minimize this function with respect to T for each point in the parameter space A = 
(n e , gioB n G), i.e. T is an empirical parameter determined by the \t minimization for each A rather than being fixed 
at the standard value To = 2.725 ± 0.002 K In principle, one should marginalize also over the galactic foreground 
spectrum |33j . However, this is a subleading effect relative to the spectral deformation caused by the photon-axion 
conversion. 

In Fig. we show the exclusion contour in the plane of n e and gioB n Q. The region above the continuous curve is the 
excluded region at 95% CL, i.e. in this region the chance probability to get larger values of xt is lower than 5%. We 
also show the corresponding 99% CL contour which is very close to the 95% contour so that another regression method 
and/or exclusion criterion would not change the results very much. Within a factor of a few, the same contours also 
hold if one varies the domain size s within a factor of 10. Comparing this exclusion plot with the iso-dimming curves 
of Fig. [3] we conclude that the entire region n e < 10~ 9 cm -3 is excluded as a leading explanation for SN dimming. 

A few comments are in order. Intergalactic magnetic fields probably are a relatively recent phenomenon in the 
cosmic history, arising only at redshifts of a few. As a first approximation we have then considered the photon-axion 
conversion as happening for present [z = 0) CMB photons. Since P 7 _* 7 is an increasing function of the photon energy 
lu, our approach leads to conservative limits. Moreover, we assumed no correlation between n e and the intergalactic 
magnetic field strength. It is however physically expected that the fields are positively correlated with the plasma 
density so that relatively high values of gioB n Q should be more likely when n e is larger. Our constraints in the region 
of n e > 10~ 10 cm -3 are thus probably tighter than what naively appears. 

V. QSO CONSTRAINTS 

CMB limits are nicely complementary to the ones obtained from the effects of photon-axion conversion on quasar 
colors and spectra [35| . One effect of photon-axion oscillations is that a dispersion is added to the quasar spectra due 
to the energy dependence of the effect. By comparing the dispersion in observed in quasar spectra with the dispersion 
in simulated ones, one can find out if the model behind each simulation is allowed. The SuperNova Observation 
Calculator (SNOC) j3|J was used j3j| to simulate the effects of photon-axion oscillations on quasar observations 
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FIG. 8: Simulated quasar spectra at z = 1 for different photon-axion oscillation scenarios. (Figure from Ref. jSSj with permission.) 
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FIG. 9: Exclusion plot for photon-axion conversion. The region to the right of the dot-dashed line is excluded by requiring achromaticity 
of SN la dimming. The region inside the dashed lines is excluded by the dispersion in QSO spectra. Assuming an intrinsic dispersion of 
5% in QSO spectra, the excluded region could be extended up to the dotted curve. The CMB argument excludes the entire region above 
the solid curve at 95% CL. (Figure from Ref. E^ .) 



(Fig. [SJ) . If the simulated dispersion is smaller than the observed, one cannot exclude the scenario since real quasars 
have an intrinsic dispersion. 

In Fig. [5] we superimpose the CMB exclusion contours with the schematic region excluded by quasars . The region 
to the right of the dot-dashed line is excluded by requiring achromaticity of SN la dimming 30] . The region inside 
the dashed lines is excluded by the dispersion in QSO spectra. Moreover, assuming an intrinsic dispersion of 5% in 
these spectra, the excluded region could be enlarged up to the dotted lines. The CMB argument excludes the region 
above the solid curve at 95% CL. 



1 We use the exclusion regions of astro-ph/0410501vl. In the published version 1.351 . corresponding to as tro- ph/0410501 j2, the iso-dimming 
curves were erroneously changed. The difference is that in version 1 the angle a in Eq. (3) of Ref. |35j[ that characterizes the random 
magnetic field direction was correctly taken in the interval 0—360° whereas in version 2 it was taken in the interval 0-90° (private 
communication by the authors). 
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A cautionary remark is in order when combining the two constraints. As we have discussed in the previous section, 
CMB limits on photon-axion conversion are model independent. On the other hand, the limits placed by the QSO 
spectra may be subject to loop holes, since they are based on a full correlation between the intergalactic electron 
density and the magnetic field strength, which is reasonable but not well established observationally. 

VI. CONSTRAINTS FROM ANGULAR DIAMETER DISTANCE 

We now turn briefly to two other types of constraint on the photon-axion conversion mechanism. The first is based 
on angular diameter distance measurements of radio-galaxies. For a source of linear radius r and angular diameter 9, 
the angular diameter distance is 

dA = y ■ (31) 

In metric theories where photons travel on null geodesies and their number is conserved, the angular distance dA and 
the luminosity distance d^ are fundamentally related by the reciprocity relation |37| 

d L {z) = (l + z) 2 d A (z) . (32) 

Photon-axion conversion in intergalactic magnetic fields would not affect the angular-diameter distance |38l l39j| and 
hence would cause a fundamental asymmetry between measurements of <iz,(z) and cIa(z)- 

In a first search for a violation of the reciprocity relation, a joint analysis of high-redshift SNe la [<i£,(z)] and radio 
galaxies [(^(z)] was undertaken [3^|. The results do not favour the loss of photons and hence disfavour mixing. 
However, this constraint is less robust than the QSO one because it is affected by possibly large systematic errors 
that are difficult to quantify |40j| . 

Since angular-diameter distance is immune to the loss of photons, the axion-conversion versus accelerating-universe 
ambiguity in the interpretation can be resolved |4l| by combining CMB acoustic peak measurements with the recent 
detection of baryon oscillations in galaxy power spectra 42] . This combination excludes a non-accelerating dark-energy 
species at the 4<r level regardless of the level of the axion coupling. 

VII. CONCLUSIONS 

We have reviewed the intriguing and phenomenologically rich mechanism of conversion of photons into very low- 
mass axion-like particles in the presence of intergalactic magnetic fields. We have examined the existing astrophysical 
and cosmological limits on this model, coming from the distortion of the CMB spectrum, from the quasar dispersion, 
and from the angular diameter distance, including the baryon oscillations detected in large-scale structure surveys. 

In particular, we have shown that the resulting CMB spectral deformation excludes a previously allowed parameter 
region corresponding to very low densities of the intergalactic medium (IGM). These limits are complementary to 
the ones derived from QSO dispersion which place serious constraints on the axion-photon conversion mechanism, 
especially for relatively large densities of the IGM. As a result, it appears that the photon-axion conversion will not 
play a leading role for the apparent SN la dimming. 

It may still happen that ultra-light or massless axions play an important cosmological role. For example, it was 
shown that by adding a photon-axion conversion mechanism on top of a dark energy model with w > — 1, one can 
mimic cosmic equations of state as negative as w ~ —1.5 |44|. Although at present there is no need for such an 
extreme equation of state, it is an interesting possibility to keep in mind, especially since alternative explanations 
as ghost /phantom fields usually pose a threat to very fundamental concepts in general relativity and quantum field 
theory. 
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APPENDIX A: PHOTON-AXION CONVERSION IN A RANDOM BACKGROUND 

In the following, we derive Eq. I|17|) along the lines of Ref. [2^ . It is assumed that photons and axions traverse N 
domains of equal length s. The component of the magnetic field perpendicular to the direction of flight Bt is constant 
within each domain and of equal strength (B = |-Bt|) m each domain, but it is assumed to have a random orientation 
in each cell. 

We begin with an initial state that is a coherent superposition of an axion and the two photon states |^4i 2) that 
correspond respectively to photons polarized parallel and perpendicular to the magnetic field in the first domain, 



Kl (0)|A 1 ) + K2 (0)|A 2 )+ Ka (0)|a) . (Al) 



The initial photon and axion fluxes are 



J 7 (0) ~ M0)| 2 + M0)| 2 , (A2) 
I a {0) ~ M0)| 2 , (A3) 

respectively. In the n-th domain the transverse magnetic field Bj~ is tilted by an angle 7„ compared to the magnetic 
field in the first domain 

|Af) = c^A,} + s n \A 2 ) , (A4) 

\Al) = - Sn |A 1 )+c ri | J 4 2 ) , (A5) 

or 

CiO) = c„/cjf(z) -SnKKz) , (A6) 

C2O) = s„«f(z) +c n K r l{z) , (A7) 

where c„ = cos7„ and s n = sin7„. Only photons polarized parallel to the magnetic field mix with axions. The values 
of the transition elements are equal in each domain since the magnitude of the magnetic field B has been assumed 
to be the same everywhere. The transition probability Pq for photon to axion oscillation in one domain is given by 
Eq. (|15|l . and the photon survival probability is 1 — Pq. At the end of the n-th domain, the photon and axion fluxes 
are 

I 7 (n + 1) ~ (l-P 4)Mz„)| 2 (A8) 

+ (1 - P S 2 n )\K 2 (z n )\ 2 + P Q \ Ka (z n )\ 2 + ... 

J„(n + 1) ~ P c 2 n \ Kl (z n )\ 2 (A9) 

+P S 2 n \K 2 ( Zn )\ 2 + (1 - P Q )\K a (z n )\ 2 + ... 

where the dots represent terms that are proportional to c n , s n , or c n s n . We have defined z n = (n — l)s. The 
coefficients K\, k 2 and taken at the beginning of the n-th domain. 

Next, we assume that the transition probability in one domain is small, i.e. Pq <C 1, and that the direction of the 
magnetic field is random, i.e. 7„ is a random variable so that 7„+i — 7„ is of order unity. Due to the randomness of 
the magnetic field, in this limit c 2 and s 2 can be replaced by their average value 1/2, while the interference terms c n , 
s n and c n s n are averaged to zero. Using 

J 7 (n) - Mz„)| 2 + | K2 (z„)| 2 , (A10) 
I a - K(z n )\ 2 , (All) 

ttti)) - i^- u i:v)(7^i (A12) 



one arrives at 





2-2(1- §P ) " 
2(l-|^o)" +1 

As the number of domains is large one can replace (1 — 3Po/2)™ +1 with the limiting function exp [— 3Poz/(2s)] to 
arrive at the final expressions 

Ly{z) = 7 7 (0)-P 7 ^ o [7 7 (0)-27 a (0)], (A13) 
I a (z) = 7 a (0)+P 7 ^ a [J 7 (0)-2/ a (0)], (A14) 
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with 



p 

1 7 — >a 



1 — exp 



3P Q z 
2s 



(A15) 
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